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ABSTRACT 



Context. The mechanisms producing the diffuse radio emission in galaxy clusters, and in particular their connection with cluster 
mergers, are still debated. 

Aims. We seek to explore the internal dynamics of the cluster Abell 1758N, which has been shown to host a radio halo and two relics, 
and is known to be a merging bimodal cluster. 

Methods. Our analysis is mainly based on new redshift data for 137 galaxies acquired at the Telescopio Nazionale Galileo, only four 
of which have redshifts previously listed in the literature. We also used photometric data from the Sloan Digital Sky Survey and from 
the Canada-France-Hawaii Telescope archive. We combined galaxy velocities and positions to select 92 cluster galaxies and analyzed 
the internal cluster dynamics. 

Results. We estimate a cluster redshift of (z> = 0.2782 and quite a high line-of-sight (LOS) velocity dispersion cry ~ 1300 km s"'. 
Our 2D analysis confirms the presence of a bimodal structure along the NW-SE direction. We add several pieces of information to the 
previous merging scenario: the two subclusters (here A1758N(NW) and A1758N(SE)) cannot be separated in the velocity analyses 
and we deduce a small LOS velocity difference (AVif.Los~ 300 km s"' in the cluster rest-frame). The velocity information successfully 
shows that A1758N is surrounded by two small groups and active galaxies infalling onto, or escaping from, the cluster. Removing the 
two groups, we estimate crv,NW ~ 1000 km s"' and cTv.se ~ 800 km s"' for A1758N(NW) and A1758N(SE), respectively. We find 
that Abell 1758N is a very massive cluster with a range of M = 2 - 3 xlO'^ h^^ Mq, depending on the adopted model. 
Conclusions. As expected for clusters that host powerful, extended, diffuse radio emissions, Abell 1758N is a major cluster merger 
just forming a massive system. 

Key words. Galaxies: clusters: individual: Abell 1758 - Galaxies: clusters: general - Galaxies: kinematics and dynamics 



1. Introduction 

Merging processes constitute an essential ingredient of the evo- 
lution of galaxy clusters (Feretti et al. 2002bi. An interesting 
aspect of these phenomena is the possible connection between 
cluster mergers and extended, diffuse radio sources: halos and 
relics. The synchrotron radio emission of these sources demon- 
strates the existence of large-scale cluster magnetic fields and 
of widespread relativistic particles. Cluster mergers have been 
proposed to provide the large amount of energy necessary for 
electron reacceleration to relativistic energies and for magnetic 
field amphfication (Tribble [19931 Feretti [T9991 Feretti I2002al 
Sarazin 2002). Radio relics ("radio gischts" as they were called 
by Kempner et al. 2004), which are polarized and elongated ra- 
dio sources located in the cluster peripheral regions, seem to be 
directly associated with merger shocks (e.g., Ensslin et al. 119981 
Roettiger et al. [T9991 Ensslin & Gopal-Krishna 120011 Hoeft et 
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al. 120041 ). Radio halos are more likely associated with the tur- 
bulence following a cluster merger, although the precise radio 
formation scenario remains unclear (re-acceleration vs. hadronic 
models, e.g., Brunetti et al. .2009 ; Ensslin et al. 2011) . Recent 
semi-analytical calculations in the framework of the turbulent 
re-acceleration scenario have allowed the community to derive 
the expectations for the statistical properties of giant radio ha- 
los, in agreement with present observations that halos are found 
in very massive clusters (Cassano & Brunetti 2005 , Cassano et 
al. 120061 ). Alternative models have been presented in the frame- 
work of hadronic models, where the time-dependence of the 
magnetic fields and of the cosmic ray distributions is taken into 
account to explain the observational properties of both halos and 
(most) relics (Keshet & Loeb .2010j . In these models the intra- 
cluster medium (ICM hereafter) magnetization is triggered by a 
merger event, in part but probably not exclusively in the wake 
of merger shocks. Unfortunately, one has been able to study 
these phenomena only recently on the basis of a sufficient statis- 
tics, i.e. a few dozen clusters hosting diffuse radio sources up 
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to z ~ 0.5 (e.g., Giovannini et al. 119991 see also Giovannini & 
Ferett i l2002l Feretti EOO? Venturi et al. '2008' Bonafede et al. 
120091 Giovannini et al. 2009). It is expected that new radio tele- 
scopes will highly increase the statistics of diffuse sources (e.g., 
LOFAR, Cassano et al.i2010ai). 

From the observational point of view, there is growing evi- 
dence of the connection between diffuse radio emission and clus- 
ter mergers, since up to now diffuse radio sources have been de- 
tected only in merging systems (see Cassano et al. I2010bl l. In 
most cases the cluster dynamical state has been derived from X- 
ray observations (Schuecker et al. 1200 1 1 Buote 120021 Cassano 
et aI. I20I0b1 l. Optical data are a powerful way to investigate the 
presence and the dynamics of cluster mergers, too (e.g., Girardi 
& Biviano l2002l l. The spatial and kinematical analysis of mem- 
ber galaxies allow us to reveal and measure the amount of sub- 
structure, and to detect and analyze possible pre-merging clumps 
or merger remnants. This optical information is indeed comple- 
mentary to the X-ray information because galaxies and the ICM 
react on different timescales during a merger (see, e.g., the nu- 
merical simulations by Roettiger et al. '1997V In this context, 
we are conducting an intensive observational and data analysis 
program to study the internal dynamics of clusters with diffuse 
radio emission by using member galaxies (DARC - Dynamical 
Analysis of Radio Clusters - project, see Girardi et al. I^O'U)- 

During our obsei-vational program, we have conducted an in- 
tensive study of the cluster Abell 1758 (hereafter A1758). A1758 
is a very rich Abell cluster having Abell richness class = 3 (Abell 
et al. 11989 '). Optically, the cluster is classified as Bautz-Morgan 
class III (Abell et al. 1989) and a Rood-Sastry irregular cluster 
"F", i.e. it has a flattened configuration (Struble & Rood 1987). 
While A1758 was classified as a single cluster by Abell, ROSAT 
images show that there are two distinct clusters (Abell 1758N 
and AbeU 1758S, hereafter A1758N and A1758S), both highly 
disturbed systems, separated by approximately 2 h^^^ Mpc along 
the N-S direction (Rizza et al. 119981 1 and both very X-ray lumi- 
nous: Lx(0. 1-2.4 keV)-i i '^« - ^'^'^ '--^ ■ 

^50 erg s 
mology) for A1758N and A1758S, respectively. 

The reference X-ray study for A1758 is that of David & 
Kempner (120041) . based on both Chandra and XMM-Newton 
data. According to their analysis, A1758N and A1758S do not 
reveal any sign of interaction between the two clusters (see also 
Durret et al. 120111 ) and their LOS recessional velocity, as mea- 
sured from the X-ray spectra, is within 2100 km s"^ suggesting 
that they likely form a gravitationally bound system. Moreover, 
David & Kempner (2004) found that A1758N is in the late stages 
of a large impact parameter merger between two hot (kTx ~1 
keV) subclusters (NW and SB subclumps), with the two remnant 
cores separated in projection by 800 /i^^q kpc and surrounded 
by hotter gas {kTx ~9-12 keV) that was probably shock-heated 
during the early stages of the merger. In particular, the Chandra 
image suggests that the X-ray NW subclump is currently moving 
toward the N, while the X-ray SB subclump is moving toward the 
SB. 

The optical luminosity distribution of the cluster red- 
sequence galaxies reveals two luminous NW and SB subclumps, 
too (see Okabe & Umetsu 2008). However, while the peak posi- 
tion of the X-ray NW subclump coincides with that of the lumi- 
nous brightest cluster member lying at NW, the X-ray SB sub- 
clump is offset by about 290 /i^,} kpc to the NW of the optical SB 
structure. The bimodal feature of A1758N in the weak-lensing 
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mass maps reported by Dahle et al. (12002) and Okabe & Umetsu 
(2008) shows that the mass and light are similarly distributed 
in A1758N. Both studies found arc candidates around the two 
mass peaks. The SB structure in the mass/galaxy map is located 
in front of the X-ray SB subclump and moves toward the south- 
east, while the NW structure shows no significant offset among 
the galaxy, ICM, and mass distributions (Okabe & Umetsu l20081 
Ragozzine et al. 120121 ). 

More recently, DuiTet et al. (1201 1) confirmed the complex 
structure of A1758N and the high core X-ray temperature (6-7 
keV) with XMM-Newton data. They also found two elongated re- 
gions of high metallicity in A1758N, likely left by ram-pressure 
stripping during the merger of two subclusters. Moreover, us- 
ing Canada-France-Hawaii Telescope (CFHT) images, Durret et 
al. (i2011 j) computed A1758N and A1758S luminosity functions 
and discussed differences from a Schechter function in the view 
of the cluster internal dynamics. 

As for the diffuse radio emission, the first evidence of a dif- 
fuse radio source in A1758N came from Kempner & Sarazin 
(120010 . Giovannini et al. (2006) first reported the clear detec- 
tion of a radio halo in the center of the cluster. More recently, 
Giovannini et al. (2009) identified a central diffuse radio emis- 
sion (halo), 0.63 /i^y Mpc in size, and two brighter peripheral 
structures on the opposite sides with respect to the cluster center 
that resemble relic radio sources (see their Fig. 7 and our Fig. [1] 
for a multiwavelength view of the cluster). A1758N is one of the 
very few clusters with a central halo and two peripheral relics 
similar to RXCJ1314.4-2515 (see Feretti et al. 2005 ). 

Spitzer/MIY'S 24 jjm observations of A1758 show a very ac- 
tive star formation in the A1758N galaxies, suggesting that dust- 
obscured activity in clusters is triggered by the details of cluster- 
cluster mergers, too (Haines et al. 12009) . 

To date, only a small amount of redshift data have been avail- 
able for A1758, they indicate a value of z ~ 0.28. No internal dy- 
namical analysis based on member galaxies has been published. 

Our new spectroscopic data for A1758N come from the 
Telescopio Nazionale Galileo (TNG). Our present analysis is 
also based on public optical photometric data from the Sloan 
Digital Sky Survey (SDSS) and the CFHT archive. 

This paper is organized as follows. We present our new opti- 
cal data and the cluster catalog in Sect. 2. We present our results 
about the cluster structure in Sect. 3. We briefly discuss our re- 
sults and present our conclusions in Sect. 4. 

Unless otherwise stated, we indicate errors at the 68% con- 
fidence level (hereafter c.l.). Throughout this paper, we use 
Ho - 70 km s ' Mpc ' in a flat cosmology with Qq - 0.3 and 
Qa = 0.7. In the adopted cosmology, 1' corresponds to ~ 254 
/z^Q kpc at the cluster redshift. 



2. Redshift data and galaxy catalog 

Multi-object spectroscopic observations of A1758N were car- 
ried out at the TNG, a 4m-class telescope, in May 2008 and May 
2009. We used DOLORBS/MOS with the LR-B Grism 1, yield- 
ing a dispersion of 187 A/mm. The detector is a 2048 x 2048 
pixels B2V CCD, with a pixel size of 13.5 //m. In total, we ob- 
served four MOS masks (one in 2008 and three in 2009) for a 
total of 146 slits. The targets were chosen within the pre-imaged 
area (a rectangle of ~8.6'xl3.7' centered on the cluster and with 
the major side clockwise rotated by 40°) after a rough color pre- 
selection. We acquired three exposures of 1800 s for each mask. 
Wavelength calibration was performed using helium -i-argon and 
helium-Hmercury-i-neon lamps. Reduction of spectroscopic data 
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Fig. 1. Multiwavelength picture of the cluster A1758N (north at the top and east to the left). Optical density peaks detected from 
our analysis of the galaxy distribution (big blue circles; see Sect. 3.4) are superimposed on the SDSS r'-band image of the cluster 
region. Labels indicate the IDs of galaxies cited in the text. Thin red contour levels show the ICM distribution we derived from 
the X-ray Chandra archival image ID 2213 (photons in the energy range 0.3-7 keV; see David & Kempner 120041 for an extensive 
analysis of these X-ray data). From Giovannini et al. (2009) we also reproduce the radio contour levels of their VLA radio image at 
1.4 GHz (thick black contours, HPBW=45"x45") and highlight with arrows the positions of the radio halo (H, in the center) and 
the two peripheral radio relics (Rl and R2). The insert on the top left is a zoom on the region of the NAT radio galaxy 1330-H507 
(ID 105), with the radio contours taken from VLA high-resolution radio images presented by O'Dea & Owen (11985) . 



was carried out using the IRAFQ package. Radial velocities 
were determined using the cross-correlation technique (Tonry & 
Davis rT979] l that is implemented in the RVSAO package (devel- 



^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



oped at the Smithsonian Astrophysical Observatory Telescope 
Data Center). Each spectrum was correlated against six tem- 
plates for a variety of galaxy spectral types: E, SO, Sa, Sb, Sc, and 
Ir (Kennicutt I 19921) . The template producing the highest value of 
•R, i.e., the parameter given by RVSAO and related to the signal- 
to-noise ratio of the correlation peak, was chosen. Moreover, all 
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Fig. 2. SDSS r'-band image of A1758N (north at the top and east to the left) showing galaxies with measured spectroscopic 
redshifts. Circles and boxes indicate cluster members and nonmember galaxies, respectively (see Table [!}. Labels indicate the IDs 
of galaxies cited in the text. 



spectra and their best correlation functions were examined visu- 
ally to verify the redshift determination. 

In eleven cases (IDs. 1, 6, 14, 31, 48, 49, 89, 110, 113, 122 
and 131; see Table [1]), we had to rely on the EMSAO package to 
obtain an estimate of the redshift. 

The formal errors as given by the cross-correlation are 
known to be smaller than the true errors (e.g., Malumuth et al. 
[1992 ; Bardelli et al. [T994l ElHngson & Yee [T994l Quintana et 
al. [2000). Duplicate observations for the same galaxy allowed 
us to estimate the true intrinsic eiTors in data of the same qual- 
ity taken with the same instiTiment (e.g. BaiTena et al. 2009J. 
Here we have double determinations for six galaxies, therefore 
we decided to apply the procedure already applied in Barrena 



et al. (2009'), obtaining that true errors are larger than formal 
cross-correlation eiTors by a factor of 2.2. For the galaxies with 
two redshift estimates, we used the mean of the two measure- 
ments and the coiTesponding eiTors. As for the radial velocities 
estimated through EMSAO we assumed the largest between the 
formal error and 100 km s ' . 

Our spectroscopic catalog lists 137 galaxies in the field of 
A1758N (see Fig.|2]i. The median eiTorin cz is 76 km s"'. 

Out of 137 galaxies, only four have a previously measured 
redshift as listed by NED. In particular, we looked at SDSS 
spectroscopic data (Data Release 7) finding eight "likely" mem- 
ber galaxies with measured spectroscopic redshift in the field 
we sampled with the TNG. Out of these, only three galaxies are 
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in common with our sample. The three redshifts agree with our 
estimates within l-2cr and the difference Vour - vsdss ranges be- 
tween 72-92 km s ' for all three galaxies. Owing to the small 
improvement we decided to not add the other five, preferring to 
work with a homogeneous sample. However, in our discussions 
we also consider the BCG of A1758S (R.A. = 13'>32™32 ?96 and 
Dec.=+50°25'02.5" -J2000.0,cz = 81802+52, r'=17.13; here- 
after BCGS). 

We used public photometric data from the SDSS (Data 
Release 7). In particular, we used g', r', and /' magnitudes, 
already corrected for the Galactic extinction, and considered 
galaxies within a radius of 25' from the cluster center. 

Table[T]lists the velocity catalog (see also Fig.|2]l: identifica- 
tion number of each galaxy and members, ID and IDm (Cols. 1 
and 2, respectively); right ascension and declination, a and 6 
(J2000, Col. 3); r' SDSS magnitude (Col. 4); heliocentric radial 
velocities, v - czq (Col. 5) with eiTors, Av (Col. 6). 

We also used public photometric data from the CFHT 
archive taken with the Megaprime/Megacam camera. We re- 
trieved from the CADC Megapipe archive (Gwyn 2009 1 the cat- 
alogs for the images in the g^'=sa ^jj^j ^Mega ^and^ and coiTected 
the magnitudes for the Galactic extinction. The total area cov- 
ered by the images was 1 .05 x 1 . 1 6 deg^ . The estimated limiting 
magnitudes (at the 5cr c.l.) are g'^'sa = 27.0 and r^'^^'' = 26.8. 

No evident unique dominant galaxy is present in A1758N, 
instead there are two BCGs: one in the NW subcluster (ID. 55, 
/ = 17.15, hereafter BCGN(NW)) and one in the SE subcluster 
(ID. 101, r' = 17.02, hereafter BCGN(SE)). They are separated 
by ~ 0.85 /i7(5 Mpc and are ~ 0.8 mags brighter than other bright 
galaxies in our sample. Moreover, Durret et al. (201 li found that 
they are characterized by a surface brightness profile flatter than 
those of other bright galaxies and suggested that the two BCGs 
are really dominant galaxies in their respective subcluster. 



2. 1 . Radio and X-ray emitting gaiaxies in the field of A1 758N 

Our spectroscopic catalog lists some radio galaxies observed in 
the sky region of A1758N. At -14" NE of BCGN(SE) we iden- 
tify the galaxy ID 105 with the radio source 1330-1-507 (O'Dea 
& Owen|T985] see our Fig.lTJ, one of the most powerful narrow- 
angle tail (NAT) radio sources known with the tail pointing to 
the SE. Rizza et al. (2003) used high-resolution VLA data to 
list a total of 11 radio sources in the field of A1758, some of 
them members of A1758N. Interestingly, source #1 of Rizza et 
al. (see flieir Table 3) is our ID 55 (the BCGN(NW)). Other ra- 
dio sources are ID 51 (#2 of Rizza et al.), ID 33 (#3, background 
galaxy), ID 124 (#6, the third-brightest galaxy in our spectro- 
scopic sample) and ID 136 (#10). 

According to Hart et al. (2009), BCGN(NW) is also a 
pointlike X-ray source. A look at the original (non-smoothed) 
Chandra X-ray image reveals that ID 124 is a faint X-ray source, 
too. 



3. Analysis of the spectroscopic sample 

3.1. Member selection 

To select cluster members among the 137 galaxies with redshifts, 
we followed a two-step procedure. We first used the ID adaptive- 
kernel method (hereafter DEDICA, Pisani 1993] and [T9961 see 



Table 1. Velocity catalog of 137 spectroscopically measured 
galaxies in the field of the cluster A1758N. f and t highlight 
the IDs of the BCGN(NW) and BCGN(SE), respectively. 



ID IDm 



tt, (5 (J2000) 



V 

(km 



Av 



1 


- 


13 32 13.95, +50 36 23.2 


21.75 


125775 


71 


2 


- 


13 32 14.70, +50 34 524 


19.65 


99042 


82 


3 


- 


13 32 15.38, +50 35 18.1 


20.01 


85168 


104 


4 


1 


13 32 16.88, +50 33 19.3 


19.79 


80993 


94 


5 


- 


13 32 17.62, +50 37 27.2 


19.11 


41801 


173 


6 


- 


13 32 17.85, +50 33 39.6 


20.10 


32788 


100 


7 


- 


13 32 17.92, +50 37 39.5 


20.18 


88373 


128 


8 


- 


13 32 18.32, +50 36 47.6 


20.23 


98217 


93 


9 


2 


13 32 18.98, +50 36 044 


21.70 


84170 


129 


10 


- 


13 32 19.28, +50 36 03.8 


20.21 


85959 


128 


11 


- 


13 32 19.45, +50 34 11.1 


20.31 


112194 


60 


12 


- 


13 32 22.27, +50 35 08.9 


19.89 


98001 


50 


13 


- 


13 32 22.42, +50 34 48.9 


19.90 


97790 


69 


14 


- 


13 32 22.47, +50 37 04.5 


21.95 


97872 


100 


15 


- 


13 32 22.59, +50 35 27.4 


20.17 


98137 


88 


16 


3 


13 32 22.89, +50 36 04.7 


20.36 


84969 


60 


17 


4 


13 32 25.27, +50 34 03.5 


18.32 


82720 


48 


18 


- 


13 32 25.98, +50 38 18.2 


20.16 


85091 


114 


19 


5 


13 32 27.59, +50 33 31.4 


20.41 


84801 


98 


20 


6 


13 32 27.59, +50 37 44.5 


19.84 


82907 


92 


21 


7 


13 32 29.23, +50 34 10.0 


19.73 


83370 


92 


22 


- 


13 32 29.57, +50 35 25.8 


20.33 


74891 


116 


23 


8 


13 32 29.93, +50 37 21.5 


20.00 


82044 


54 


24 


- 


13 32 30.23, +50 36 37.1 


20.32 


113677 


79 


25 


9 


13 32 30.27, +50 33 49.0 


19.89 


82680 


70 


26 


10 


13 32 30.62, +50 36 25.5 


20.04 


82111 


66 


27 


- 


13 32 31.64, +50 36 17.7 


19.73 


98544 


76 


28 


11 


13 32 31.88, +50 34 49.9 


20.37 


83238 


123 


29 


- 


13 32 32.64, +50 34 27.9 


18.55 


53251 


53 


30 


12 


13 32 32.86, +50 36 04.5 


19.96 


85335 


62 


31 


- 


13 32 33.31, +50 34 02.5 


21.70 


187664 


100 


32 


- 


13 32 33.72, +50 31 03.4 


19.90 


112311 


98 


33 


- 


13 32 33.86, +50 30 50.0 


19.01 


112310 


75 


34 


- 


13 32 34.12, +50 30 39.7 


20.22 


113091 


116 


35 


13 


13 32 34.31, +50 33 03.1 


19.47 


87495 


50 


36 


14 


13 32 34.41, +50 32 11.3 


19.78 


80234 


41 


37 


15 


13 32 34.49, +50 33 18.5 


19.11 


84766 


70 


38 


16 


13 32 34.57, +50 33 16.2 


21.48 


83263 


223 


39 


17 


13 32 34.69, +50 31 39.1 


20.45 


84767 


92 


40 


18 


13 32 34.70, +50 33 24.7 


19.43 


81303 


72 


41 


19 


13 32 34.75, +50 31 43.3 


20.46 


86327 


79 


42 


20 


13 32 34.96, +50 32 37.3 


18.75 


84194 


44 


43 


- 


13 32 35.07, +50 38 30.5 


19.68 


98294 


119 


44 


21 


13 32 35.21, +50 32 36.2 


18.25 


81900 


71 


45 


22 


13 32 35.69, +50 33 38.1 


20.29 


83293 


54 



^ For a comparison between Megacam 

and SDSS filters, check the web page 

|http://www2.cadc-ccda.hia-iha.nrc-ciirc.gc.ca/megapipe/docs/filters.html 



also Fadda et al. 1 19961 Girardi et al. 1 19961 ). We searched for sig- 
nificant peaks in the velocity distribution at >99% c.l.. This pro- 
cedure detects A1758N as a peak at z ~ 0.278 populated by 104 
galaxies considered as candidate cluster members (in the range 
79 373 < V < 88 373 km s"', seeFig.[3]l. Out of 33 nonmembers, 
8/25 are foreground/background galaxies. 

All galaxies assigned to the cluster peak were analyzed in 
the second step, which combines position and velocity informa- 
tion, i.e., the "shifting gapper" method by Fadda et al. ( 119961 ). 
This procedure rejects galaxies that are too far in velocity from 
the main body of galaxies within a fixed bin that shifts along the 
distance from the cluster center. The procedure is iterated un- 
til the number of cluster members converges on a stable value. 



W. Boschin et al.: Abell 1758N: new optical insights on a merging cluster with diffuse radio emission 
Table 1. Continued. Table 1. Continued. 



ID IDm 



Q',5(J2000) 



V 

(km 



Av 



46 


- 


13 32 35.73 


+50 34 49.8 


20.83 


113693 


82 


47 


23 


13 32 35.84 


+50 33 52.0 


20.50 


83709 


76 


48 


- 


13 32 36.46 


+50 35 50.6 


20.43 


87368 


100 


49 


- 


13 32 36.88 


+50 30 37.7 


20.50 


112681 


100 


50 


24 


13 32 36.89 


+50 32 59.7 


19.54 


83963 


63 


51 


25 


13 32 37.58 


+50 33 05.7 


19.35 


82279 


44 


52 


26 


13 32 38.04 


+50 30 57.2 


20.58 


84344 


79 


53 


27 


13 32 38.33 


+50 31 30.3 


20.92 


86862 


255 


54 


28 


13 32 38.33 


+50 33 30.1 


20.48 


84176 


120 


55 


t29 


13 32 38.41 


+50 33 35.7 


17.15 


83489 


57 


56 


30 


13 32 38.46 


+50 3141.4 


20.59 


85204 


113 


57 


31 


13 32 38.48 


+50 32 53.0 


20.14 


83376 


76 


58 


32 


13 32 38.58 


+50 33 43.9 


19.52 


83791 


79 


59 


33 


13 32 39.43 


+50 34 45.0 


18.97 


83208 


50 


60 


34 


13 32 39.54 


+50 33 45.0 


20.13 


83880 


76 


61 


35 


13 32 39.54 


+50 34 32.0 


19.01 


87769 


85 


62 


36 


13 32 39.76 


+50 33 14.4 


20.40 


82872 


128 


63 


37 


13 32 39.80 


+50 32 41.0 


19.37 


83803 


84 


64 


38 


13 32 39.94 


+50 32 23.9 


20.94 


83497 


172 


65 


39 


13 32 40.35 


+50 34 49.4 


19.50 


85818 


70 


66 


40 


13 32 40.51 


+50 35 39.6 


18.56 


81400 


47 


67 


41 


13 32 40.52 


+50 32 16.0 


21.97 


85932 


182 


68 


42 


13 32 40.59 


+50 33 15.4 


19.63 


80255 


54 


69 


43 


13 32 40.96 


+50 33 46.3 


17.93 


83667 


38 


70 


44 


13 32 41.47 


+50 32 10.7 


20.38 


80544 


57 


71 


45 


13 32 42.03 


+50 34 34.7 


20.19 


87429 


132 


72 


46 


13 32 42.21 


+50 32 30.1 


20.56 


82939 


85 


73 


47 


13 32 42.35 


+50 32 33.1 


20.94 


81072 


91 


74 


48 


13 32 43.30 


+50 29 26.2 


18.99 


82769 


63 


75 


49 


13 32 43.43 


+50 33 28.8 


18.80 


85478 


60 


76 


50 


13 32 43.77 


+50 31 14.6 


19.35 


83504 


50 


77 


51 


13 32 44.77 


+50 33 26.0 


19.78 


84664 


97 


78 


52 


13 32 44.96 


+50 31 57.2 


20.03 


84337 


82 


79 


53 


13 32 45.35 


+50 33 24.6 


19.47 


85798 


82 


80 


54 


13 32 45.94 


+50 32 04.7 


19.50 


79995 


76 


81 


55 


13 32 46.36 


+50 36 41.5 


18.88 


86035 


47 


82 


56 


13 32 46.90 


+50 3148.5 


20.48 


83825 


113 


83 


57 


13 32 47.00 


+50 32 02.0 


19.77 


84249 


76 


84 


58 


13 32 47.92 


+50 32 09.6 


20.03 


82909 


82 


85 


- 


13 32 48.47 


+50 28 06.7 


19.29 


109410 


85 


86 


59 


13 32 48.71 


+50 31 54.9 


19.94 


82757 


63 


87 


60 


13 32 48.72 


+50 31 21.7 


19.47 


83562 


63 


88 


61 


13 32 49.25 


+5031 31.1 


19.89 


82110 


88 


89 


- 


13 32 49.28 


+50 28 50.3 


20.45 


161829 


100 


90 


62 


13 32 49.32 


+50 33 55.9 


19.95 


79788 


76 



Following Fadda et al. (1 19961 1. we used a gap of 1000 km s ' - 
in the cluster rest-frame - and a bin of 0.6 /i^,j Mpc, or large 
enough to include 15 galaxies. For the center of A1758N we 
adopted the position of the BCGN(NW) [R.A. = 13h32™38 ?41, 
Dec.=+50°33'35.7" (J2000.0)].We selected this center because, 
contrary to the SB subcluster, the peak position of the X-ray NW 
subclump coincides with that of the BCGN(NW). Moreover, the 
NW subcluster is more luminous, massive, and hot than the SE 
subcluster (Okabe et al. 120081 Durret et al. 1201 H . which sug- 
gests that it is the primary component of A1758N. The "shifting 
gapper" procedure rejected another 12 interlopers. Most of them 
are not very far from the main body, but we are confident enough 
in our rejection because 6 out of 12 are emission line galaxies 
(ELGs hereafter), which increases their probability of belonging 
to the field or to the infalling cluster region (see Sect. l4.1l for a 





ID 


IDm 


tt, 5 (J2000) 


r' 


V 

(km 


Av 


91 


_ 


13 32 49.73, +50 28 41.6 


18.64 


55954 


60 


92 


63 


13 32 50.00 


+50 32 26.0 


19.98 


84288 


79 


93 


64 


13 32 50.03 


+50 33 53.8 


20.40 


81273 


60 


94 


- 


13 32 50.62 


+50 27 52.3 


19.49 


80403 


72 


95 


65 


13 32 50.72 


+50 27 46.7 


19.55 


82691 


79 


96 


66 


13 32 51.01 


+50 33 08.8 


19.89 


84526 


79 


97 


67 


13 3251.10 


+50 27 45.8 


19.70 


81557 


82 


98 


68 


13 32 51.42 


+50 33 04.7 


17.95 


85536 


57 


99 


- 


13 32 51.50 


+50 33 10.3 


20.16 


130803 


148 


100 


69 


13 32 51.79 


+50 27 13.5 


18.01 


83197 


54 


101 


m 


13 32 52.10 


+50 31 34.1 


17.02 


83864 


70 


102 


71 


13 32 52.21 


+50 31 21.8 


20.36 


83139 


72 


103 


- 


13 32 52.33 


+50 27 42.8 


20.60 


109756 


204 


104 


72 


13 32 52.41 


+50 31 12.8 


20.23 


82717 


75 


105 


- 


13 32 52.93 


+50 3146.1 


18.52 


79682 


41 


106 


73 


13 32 53.10 


+50 3135.5 


19.64 


83648 


63 


107 


- 


13 32 53.24 


+50 30 29.0 


18.69 


99000 


47 


108 


74 


13 32 53.57 


+50 27 36.0 


20.42 


83346 


94 


109 


- 


13 32 53.69 


+50 30 53.3 


19.06 


98877 


104 


110 


- 


13 32 53.75 


+50 26 56.2 


20.81 


212245 


100 


111 


75 


13 32 53.92 


+50 29 57.7 


19.90 


82041 


66 


112 


76 


13 32 53.94 


+50 32 20.7 


19.13 


79836 


47 


113 


- 


13 32 54.86 


+50 30 27.8 


20.89 


52848 


100 


114 


77 


13 32 54.94 


+50 33 14.7 


19.82 


79373 


54 


115 


78 


13 32 55.00 


+50 31 27.6 


21.38 


84114 


154 


116 


79 


13 32 55.15 


+50 31 25.4 


18.85 


85175 


72 


117 


80 


13 32 55.61 


+50 31 28.5 


19.01 


85239 


47 


118 


- 


13 32 55.64 


+50 30 55.8 


18.83 


55708 


98 


119 


- 


13 32 56.08 


+50 30 17.5 


18.98 


31474 


70 


120 


81 


13 32 59.32 


+50 30 19.9 


20.66 


83353 


230 


121 


82 


13 32 59.59 


+50 29 27.6 


19.24 


82055 


50 


122 


- 


13 32 59.75 


+50 30 05.4 


20.61 


86250 


100 


123 


83 


13 33 02.05 


+50 27 17.1 


19.84 


82098 


54 


124 


84 


13 33 02.07 


+50 29 28.0 


17.92 


83651 


47 


125 


85 


13 33 02.12 


+50 30 09.2 


20.88 


83434 


76 


126 


- 


13 33 02.71 


+50 3142.3 


18.87 


80377 


82 


127 


- 


13 33 03.67 


+50 28 25.3 


19.87 


93943 


91 


128 


86 


13 33 06.72 


+50 31 24.1 


18.62 


85598 


47 


129 


- 


13 33 07.43 


+50 29 19.6 


20.35 


82058 


62 


130 


87 


13 33 07.48 


+50 32 00.7 


18.93 


82508 


47 


131 


- 


13 33 09.83 


+50 28 33.2 


19.57 


87007 


100 


132 


88 


13 33 09.90 


+50 29 02.3 


19.94 


83139 


76 


133 


89 


13 33 10.19 


+50 28 51.3 


18.92 


82496 


72 


134 


- 


13 33 10.72 


+50 31 15.6 


18.87 


86227 


40 


135 


90 


13 33 10.93 


+50 28 53.7 


19.90 


82857 


91 


136 


91 


13 33 12.99 


+50 29 01.3 


18.52 


82798 


41 


137 


92 


13 33 13.84 


+50 29 55.3 


20.77 


83782 


94 





more detailed discussion on these ELGs). Finally, we obtained a 
sample of 92 fiducial cluster members (see Fig.lUi. The removed 
interlopers are shown in Fig.|5](top panel). 



3.2. General cluster properties 

By applying the biweight estimator to the 92 cluster members 
(Beers et al. 11990 ROSTAT software), we computed a mean 
cluster redshift of (z) = 0.2782+0.0005, i.e. (v) = (83390+139) 
km s"'. We estimated the LOS velocity dispersion, cry, by us- 
ing the biweight estimator and applying the cosmological cor- 
rection and the standard correction for velocity errors (Danese et 
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Fig. 3. Redshift galaxy distribution. The solid line histogram 
refers to the 104 galaxies assigned to A1758N according to the 
DEDICA reconstruction method. 



al. 119801 1. We obtained cry - 1329^[j^ km s"', where errors are 
estimated through a bootstrap technique. 

To evaluate the robustness of the cry estimate, we analyzed 
the velocity dispersion profile (Fig.|5]l. The integral profile flat- 
tens in the external cluster regions, as found for most nearby 
clusters (e.g., Fadda et al. 1996; Girardi et al. 11996) . The sharp 
rising of the integral and differential profiles in the internal re- 
gion is discussed in Sect. 14. II 



3.3. Velocity distribution 

We analyzed the velocity distribution to search for possible devi- 
ations from Gaussianity that might provide important signatures 
of complex dynamics. For the following tests, the null hypothe- 
sis is that the velocity distribution is a single Gaussian. 

We estimated three shape estimators, i.e., the kurtosis, the 
skewness, and the scaled tail index STl. We found ST1= 1.294, 
i.e. that the velocity distribution departs from Gaussianity at the 
95%-99% c.l. (according to Table 1 of Bird & Beers 1199^ . The 
velocity distribution thus shows evidence for a heavy tailed dis- 
tribution. 

We then investigated the gaps in the velocity distribution. 
We followed the weighted gap analysis presented by Beers et 
al. ( fT99n[1992t ROSTAT software). We looked for normalized 
gaps larger than 2.25, since in random draws of a Gaussian dis- 
tribution they arise at most in about 3% of the cases, independent 
of the sample size (Wainer and Schacht fl 9781 1. We detected three 
significant gaps (at the 97%, 97%, and 98.6% c.ls), which divide 
the cluster into four groups of 13, 9, 14, and 56 galaxies from 
low to high velocities (hereafter GVl, GV2, GV3, and GV4, see 
Fig.m. Both BCGN(NW) and BCGN(SE) were assigned to the 
GV4 peak. 

Following Ashman et al. (I1994I I. we also applied the Kaye's 
mixture model (KMM) algorithm. We found that a partition of 
13, 7, 41, and 31 galaxies is a significantly more accurate de- 



Fig. 4. The 92 galaxies assigned to the cluster A1758N. Upper 
panel: Galaxy velocity distribution. The arrows indicate the ve- 
locities of the BCGN(NW) and BCGN(SE). We also show the 
velocity of BCGS of A1758S as taken from the SDSS. Lower 
panel: Stripe density plot where the arrows indicate the positions 
of the significant gaps. 



scription of the galaxy distribution than a single Gaussian (at the 
~ 97% c.l.). 



3.4. Analysis of the 2D galaxy distribution 

By applying the 2D adaptive-kernel method (2D-DED1CA) to 
the positions of the 92 member galaxies we found that the clus- 
ter is elongated in the SE-NW direction. In particular, we de- 
tected three peaks with high statistical significance: the main 
peak, which is coincident with the BCGN(NW), a secondary 
peak coincident with the BCGN(SE), and a third one, which is 
intermediate (Fig.|6]l. 

Our spectroscopic data do not cover the entire cluster field 
and are affected by magnitude incompleteness. More precisely, 
we found the spectroscopic catalog has a ~ 50% completeness 
for r' <19 and the completeness level quickly drops to <30% 
for r' >21. To overcome these problems, we used the SDSS and 
CFHT photometric data samples, which cover a larger spatial 
region. The SDSS sample has the advantage to have photometry 
available in several magnitude bands, while the CFHT sample 
allows us to extend our analysis to fainter galaxies. 

From the CFHT photometric catalog we selected likely clus- 
ter members on the basis of the (g^'^^^-p^'^s^) vs. r'^'^s^ color- 
magnitude relation (hereafter CMR), which indicate the locus of 
the red-sequence galaxies. To determine the CMR we applied the 
2cr-clipping fitting procedure to the spectroscopic cluster mem- 
bers. We obtained ^'^'=sa_^"ega^2. 121-0.041xr'^'=sa (,„ 55 galax- 
ies of our spectroscopic sample, which agrees very well with the 
CMR directly fitted on the photometric catalog recovered from 
CFHT images by DuiTet et al. d^OTTT l. Out of the CFHT photo- 
metric catalog, we considered as likely "red" cluster members 
the objects lying within 0.15 mag of the CMR, where ~ 0.15 is 
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Fig. 5. Top panel: rest-frame velocity vs. projected distance 
from the cluster center (here the BCGN(NW)). Small/black and 
large/red squares indicate the member galaxies and those re- 
jected as interlopers by the shifting gapper procedure, respec- 
tively. Red crosses indicate ELGs. The two green circles indi- 
cate the BCGN(NW) and the BCGN(SE). The large/blue circle 
indicates the BCG of A1758S using the SDSS redshift. The cyan 
triangles and rotated triangles indicate galaxies detected as mem- 
bers of HT2 and HT12 subclusters detected through the Htree- 
method (see Sect. 13.51 1. Middle and bottom panels: differential 
(big circles) and integral (small points) profiles of mean veloc- 
ity and LOS velocity dispersion, respectively. For the differential 
profiles, we plot the values for six annuli from the center of the 
cluster, each containing 15 galaxies (large green symbols). For 
the integral profiles, the mean and dispersion at a given (pro- 
jected) radius from the cluster-center is estimated by considering 
all galaxies within that radius - the first value computed on the 
five galaxies closest to the center. The error bands at the 68% c.l. 
are also shown. In the lower panel, the horizontal line represents 
the X-ray temperature (7 keV) estimated from David & Kempner 
(I2004I I for the core of the two subclusters, see also DuiTet et al. 
( I2011I I. transformed to cry assuming the density-energy equipar- 
tition between ICM and galaxies, i.e. ySspec = 1 (see Sect. |4]i. 



the error associated to the fitted intercept. Figure shows that 
our selection criterion seems adequate to select red-sequence 
galaxies, which are good tracers of the cluster substiTicture (e.g., 
Lubin et al. 2000) and, above all, avoiding nonmember galax- 
ies that might bias our 2D analysis. As for the sample with 
^Mega ^ 21^ we can check our selection of likely members: we 
can recover 75% of spectroscopic members (i.e. 68 out of 90), 
but 30% of nonmembers are also selected (13 out of 42). 

Figure[8]shows the galaxy density contours obtained with the 
2D-DEDICA method for two samples: the CFHT likely mem- 
bers with r^'^s" < 21 and those with 21 < r^'^^" < 23 (1657 
and 3288 galaxies in the whole sample; 293 and 374 galaxies 
within 2 h^^ Mpc, respectively). As for the luminous galax- 
ies, we found two very significant peaks along the SE-NW di- 
rection, corresponding to the BCGN(SE) and BCGN(NW): the 




2 0-2 

X [arcmin] 

Fig. 6. Spatial distribution on the sky and relative isodensity con- 
tour map of the 92 spectroscopic cluster members, obtained with 
the 2D-DEDICA method. The BCGN(NW) is taken as the clus- 
ter center. Crosses indicate the location of the BCGN(NW) and 
BCGN(SE). 



NW-peak is the richest. These two peaks are separated by ~0.75 
/z^o Mpc, with the SE-peak well coincident with the BCGN(SE) 
and the NW-peak at a distance of ~0.15 /i^,5 Mpc from the 
BCGN(NW). A third very significant peak lies in the South, co- 
incident with the position of the cluster A1758S. As for the re- 
sults with r' < 21, Table|2]lists the number of assigned members, 
A^s (Col. 2); the peak position (Col. 3); the density (in arbitrary 
units relative to the densest peak fixed to be 1), ps (Col. 4); the 
value of ;t^^ for each peak,;i'j (Col. 5). The second sample (lower 
panel of Fig. [HJ contains galaxies fainter by two mags: as dis- 
cussed at the end of Sect. 14.1! these galaxies seem to trace the 
whole A1758N system and not its substructure. 

When more than two colors are available, it is very effec- 
tive to select cluster galaxies in the color-color space (Goto et al. 
I2OO2I 1 and both SDSS r'-i' and g'-r' colors are generally used at 
z < 0.4 (Lu et al. 2009; Lopes 2007: his Eq. 1 and refs. therein). 
In particular, the difference between passive and star-forming 
galaxies is larger in g'-r' color than in r'-i' color (see Lu et 
al. 120091 their Fig. 6) but r'-i' color has a much smaller scatter 
(Lopes 120071 his Fig. 1). Out of the SDSS photometric catalog, 
we considered as likely "red" cluster members the objects ly- 
ing within 0.1 and 0.2 mag of the CMRs (/-;' vs. r') and {g'-r' 
vs. r'), respectively. When comparing this member selection to 
that performed on the CFHT catalog, we obtain a comparable re- 
sult with a modest, better rejection of nonmembers. Indeed, we 
recover 67 out of 88 spectroscopic members and select only 10 
out of 42 nonmembers (when using the g'-r' cut alone we obtain 
67 out of 88 and 14 out of 42). As for the galaxy density map, the 
SDSS result confirms the result obtained with the bright CFHT 
sample: the NW-peak is still richer than the SE-peak (77 vs. 52 
galaxies), although the SE-peak is now the densest (by a factor 
1.5). 



W. Boschin et al.: Abell 1758N: new optical insights on a merging cluster with diffuse radio emission 



1.5 



I 



t>0 



0.5 - 





, 1 , , 


, 1 , , 




' - 


= --_ 


--- — 


s % m 




- 




■'- ^^?-y"~ 


— 


13 


- 


13 


01313 


" 


- 




13 = 


13 


- 


- 




HI ° |x] 




- 


— 


m 






— 


- 




13 




- 


1 


1 


1 




- 



18 



20 



22 



pMega 



Fig. 7. CFHT gMega_^Mega ^^ ^ega diagram for galaxies with 
available spectroscopy. Small/black squares indicate member 
galaxies and large/red squares indicate nonmember galaxies. 
Red crosses indicate ELGs. The solid line gives the CMR de- 
termined on member galaxies; the dashed line is drawn at ±0. 15 
mag from this value. The green circles indicate the BCGN(NW) 
and BCGN(SE) of A1758N: these are the only galaxies of the 
sample that likely suffer for photon count saturation in the CFHT 
catalog (because they are brighter than r'^'^s^ ~ 17). 

Table 2. 2D substructure from the CFHT photometric sample. 



Subclump 



N/ Q'(J2000), (5(J2000) ps xl 
h : m : s, ° : ' : " 



2D - N(NW) (CFHT r^=s^' < 21) 
2D - S (CFHT r^-^s" < 21) 
2D - N(SE) (CFHT r^^^^" < 21) 



87 13 32 37.7 + 50 33 19 1.00 54 
69 13 32 22.5 + 50 24 39 0.25 21 
53 13 32 52.1+50 3133 0.90 51 



Values are not background-subtracted. 



3.5. 3D-analysis 

The existence of correlations between positions and velocities of 
cluster galaxies is a characteristic of true substructures . Here we 
used several approaches to perform a 3D-analysis of the cluster. 

We found no evidence for a significant velocity gradient (see, 
e.g., den Hartog & Katgert 1996i and Girardi et al. il996i for the 
details of the method). 

We computed the A-statistics devised by Dressier & 
Schectman (1988. hereafter DS-test), which is recommended by 
Pinkney et al. ( 119961 ) as the most sensitive 3D test. For each 
galaxy, the deviation 5 is defined as 6j - [(Nnn + 1)/Cv]['^^i ~ 
V)^ + ((Tv,i - cTy)^], where the subscript "1" denotes the local 
quantities computed over the A^nn = 10 neighbors of the galaxy. 
A is the sum of the 6 of the individual A^ galaxies and gives the 
cumulative deviation of the local kinematical parameters (mean 
velocity and velocity dispersion) from the global cluster param- 
eters. The significance of A, i.e. of substructure, is checked by 




5 

X [arcmin] 




Fig. 8. Spatial distribution on the sky and relative isodensity con- 
tour map of CFHT photometric cluster members with r^'=g^ < 2 1 
(upper panel) and 21 < r'^'=g^ < 23 [lower panel), obtained 
with the 2D-DEDICA method. The BCGN(NW) is taken as 
the cluster center. Northern crosses indicate the location of the 
BCGN(NW) and BCGN(SE) of A1758N. The southern cross in- 
dicates the BCGS of A1758S. For A1758N, the two large "O" 
labels are centered on the two main peaks of the mass distribu- 
tion (Okabe et al. ,2008) . 



running 1000 Monte Carlo simulations, randomly shuffling the 
galaxy velocities. We found a significant presence of substruc- 
ture at the 97% c.l.. 

Following Pinkney et al. ( 1996; see also Ferrari et al. 2003}, 
we applied two more classical 3D tests: the e-test (Bird, 1993T l 
based on the projected mass estimator and the centroid shift or 
ff-test (West & Bothun 1990). The details of these tests can be 
found in the above papers. We merely point out that we consid- 
ered ten as the number of the nearest neighbors for each galaxy 
and used Monte Carlo simulations to compute the substructure 
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Fig. 9. Spatial distribution of the 92 cluster members, each 
marked by a circle. The BCGN(NW) is taken as the cluster cen- 
ter. Upper panel: the cluster velocity field: the larger the cir- 
cle, the larger the galaxy velocity. Lower panel: the result of 
the (modified) DS-test: the larger the circle, the larger the devia- 
tion (5s,, of the local velocity dispersion from the global velocity 
dispersion. Thin/blue and thick/red circles show where the local 
velocity dispersion is lower or higher than the global value. The 
two large green squares indicate the positions of the two BCGs. 



significance. The application of the e-test leads to a significant 
presence of substructure at the 99.4% c.l.. 

To better understand the results of the DS-test described 
above, we also considered two kinematical estimators alternative 
to the 6 parameter of the DS-test, i.e. we considered separately 
the contributes of the local mean 5\ - [(M,n + l)/o"y](Vi - V)^], 

and dispersion 5s - [(A^nn + l)/o'v](^''"v,i~'^v)^] (see, e.g. Girardi 
et al. |1997 ' FeiTaii et al. '200Tl. When considering the 6s estima- 
tor, we found evidence for a peculiar local velocity dispersion 
at the 99.2% c.l.. Figure |9]- lower panel - shows the distri- 
bution on the sky of all galaxies, each marked by a circle: the 
larger the circle, the larger the deviation d^j of the local veloc- 
ity dispersion from the global cluster value. Figure |9] shows: i) 
as A1758N(NW) is well detected as a region of low cry,!, ii) a 
few galaxies in the NE central region have high values of cryj. 
The first result is expected for/inteipreted as a system having a 
relaxed core owing to circular orbits or galaxy merger phenom- 
ena (i.e., Girardi et al. [T996l and |1997 Menci & Fusco-Femiano 
11996 1). The second result deserves more attention: we compared 
the distribution of individual (J, , values of real galaxies with 
those of simulated clusters (Biviano et al. '2002V Figure [10] sug- 
gests that significantly high S^j values for real data are those with 
(5s,i > 1.8. A part from several galaxies with peculiarly low cry j, 
the 1.8 value only selects four galaxies with peculiarly high cry i 
(the largest thick/red circles in Fig. |9]- lower panel, IDs 81, 90, 
93, and 96), two of which are ELGs. 

Then we searched for a possible physical meaning of the four 
subclusters determined by the three weighted gaps. We com- 
pared the spatial galaxy distributions of GVl, GV2, GV3, and 



Fig. 10. Distribution of d^j deviations of the Dressler-Schectman 
analysis for the 92 member galaxies. The solid line represents 
the observations with Poissonian eiTors. The dashed line is the 
distribution for the galaxies of simulated clusters, normalized to 
the observed number. 



GV4 two by two. GVl galaxies are distributed in a round region 
between the NW and the SE clumps, and do not follow the elon- 
gated shape of the sample. We verified that GVl differs from 
GV3 (and from GV2-HGV3) at the 94% c.l. (95% c.l.) according 
to the 2D Kolmogorov-Smimov test (2DKS-test, Fasano et al. 
119871 see Fig. [TlTi. 

We finally resorted to the method devised by Serna & Gerbal 
(119961 hereafter Htree-method; see, e.g., Durret et al. 120 101 for 
a recent application). This method uses a hierarchical cluster- 
ing analysis to determine the relationship between galaxies ac- 
cording to their relative binding energies. The method assumes a 
constant value for the mass-to-light ratio of galaxies and Serna & 
Gerbal (1996) suggested a value comparable to that of clusters. 
Here we took a value of MjLr-lSQ hjo Mq/Lq, as suggested by 
large statistical studies (e.g., Girardi et al. i2000i Popesso et al. 
2005). 

Figure [12] shows the resulting dendogram, where the to- 
tal energy appears horizontally. Galaxy pairs and subgroups of 
galaxies appear with a lower total energy. We note that there 
are no important subsystems. In particular, the BCGN(SE) and 
BCGN(NW) he at the bottom of the potential well, although in 
two separate subsystems (a pair and a quadruplet, respectively). 
Instead, there are several small subsystems. We highlight here 
two that are situated at the highest energy levels, which might 
largely bias the computation of the global velocity dispersion. 
In Figure [12] we indicate HT2, a triplet of close galaxies at high 
velocity (v > 87 000 km s"'), and HT12, which is formed by 
seven galaxies (v ~ 80 000 km s"' ) that lie in the central cluster 
region. Note that all seven galaxies of HT12 are also members 
of the GVl group detected in our ID analysis. Other significant 
subsystems of HTll, here interpreted as the main system, are 
small and generally formed by close galaxies. We suspect that 
their detection might be very sensitive to our incomplete spa- 
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Fig. 11. Spatial distribution on the sky of the cluster galaxies, 
showing the four groups recovered by the weighted gap analysis. 
Solid/blue circles, open circles, open squares, and dots indicate 
the galaxies of GVl, GV2, GV3, and GV4, respectively. The 
BCGN(NW) is taken as the cluster center. Very large open green 
squares indicate the positions of BCGN(NW) and BCGN(SE). 
Large/red crosses indicate ELGs. 

tial sampling and we do not discuss them. Figure [T3] shows the 
spatial distributions of HT2 and HT12 with respect to the main 
system. 

4. Discussion and conclusions 

The high value of the velocity dispersion cry = 1329^|jg 
km s"' is comparable to the values found for hot, massive clus- 
ters (e.g., Mushotzky & Scharf |T997l Girardi & Mezzetti BOOTT) . 
However, our estimate of cry is somewhat higher than the X- 
ray temperature of A1758N subclusters under the assumption of 
the equipartition of energy density between ICM and galaxies. 
Indeed, assuming that kTx ~ 7 keV, we obtained ySspec ~ 1 -24 
to be compared with jSspec = 1, where jSspec = cr\l{kTxlnm^) 
with n - 0.58 the mean molecular weight and nij, the proton 
mass (see also Fig. |5]). Indeed, our cry would instead predict a 
value of kTx-lOJ keV, which would be more in line with the 
temperature found in the region beyond the two cores, where the 
gas is already shocked by the merger and the X-ray temperature 
likely already enhanced to the value of the whole forming clus- 
ter (David & Kempner 2004 and discussion therein). Below we 
discuss the bimodal nature of A1758N. 

4.1. Cluster structure and mass 

In agreement with previous results (see Sect. [T), our 2D analysis 
highlights two subclusters in the galaxy distribution about 3'(~ 
0.75 /;^Q Mpc) away from each other along the SE-NW direction, 
hereafter A1758N(NW) and A1758N(SE). 

We cannot separate A1758N(NW) and A1758N(SE) from 
the velocity information. Indeed, the small difference between 




Fig. 12. Dendogram obtained through the Sema & Gerbal (11996) 
algorithm. The abscissa is the binding energy (here in arbitrary 
units with the deepest negative energy levels on the left) while 
the catalog numbers of the various member galaxies are shown 
along the ordinate (IDm in Table [T]i. 



the LOS velocities of the two BCGs (< 300 km s ' in the rest- 
frame) is by itself an indication that the two subclusters have low 
LOS relative velocity. In this study we were able to fully exploit 
our sample of 92 cluster members. From i) the inspection of the 
galaxy velocity field in Fig. |9] (upper panel), ii) the absence of 
a velocity gradient, iii) the absence of local velocity peculiari- 
ties via the DS-test, we confirm that the velocity difference be- 
tween the two subclusters is very small. This is clearly shown by 
Fig. [14]- upper panel, where the integral mean velocities com- 
puted around A1758N(NW) and A1758N(SE) are very similar, 
lower than the velocity difference between the two BCGs. 

Taking into account the similar velocities of the two BCGs, 
the sharp increase of the velocity dispersion profile within 0.5 
/z^Q Mpc (see Fig. Island Fig.fT4l- faint red line - for both the sub- 
clusters) is not directly explainable with the presence of two sys- 
tems at different mean velocity. Rather, this might be due to the 
presence of subclusters or individual galaxies infalling onto, or 
escaping from, the cluster. Indeed, beyond the usual infall onto 
the cluster from the large-scale structure, the recent merger of 
two subclumps may result in outflying galaxies as shown in sim- 
ulations (e.g., Czoske et al . 120021 1 and seen in a few clusters (see 
e.g., the plume of outflying galaxies in Abell 3266 by Quintana 
et al. 1999 and Flores et al. 2000, the structure of C10024+1654 
by Czoske et al. l2002l l. In particular, the simulations reproducing 
C10024-H1654 by Czoske et al. (|20U2] | show that during a head- 
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Fig. 13. Spatial distribution on the sky of the cluster galaxies 
showing the subsystems recovered by the Htree analysis. Solid 
blue and red circles indicate the galaxies of HT2 and HT12, re- 
spectively. Triangles indicate the galaxies of the main system 
(HTl). The BCGN(NW) is taken as the cluster center. Very 
large/green squares indicate the positions of BCGN(NW) and 
BCGN(SE). Large/red crosses indicate ELGs. Large red dia- 
monds indicate dusty star-forming galaxies in the sample of 
Haines et al. (i2009b . 



head collision at 3000 km s"', the outer regions of the smaller 
cluster have become unbound and are streaming radially away 
from the impact location with velocities perpendicular to the en- 
counter on the order of 1000 km s ' . 

The velocity information is useful for detecting the galaxy 
groups connected with A1758N. Through our 3D Htree analysis 
we detected the high velocity HT2 group and the low velocity 
HT12 group. The HT2 and HT12 groups lie in the tails of the 
velocity distribution (which, in fact, is a heavy tailed distribu- 
tion; see Sect. l3.3l l. In particular, HT12 is a subsample of GVl, 
the only group detected in the ID analysis that also has a pecu- 
liar spatial distribution. Spatially, HT2 and HTl 1 lie in the cen- 
tral region roughly between the two BCGs, where, in the NE, 
the DS analysis detects a region of high local velocity disper- 
sion. Figure [14] - lower panel - shows the resulting profile of 
the velocity dispersion when only the main system HTl 1 is con- 
sidered: the effect is that the value of cry decreases to ~ 1000 
km s"'. 

As for individual/peculiar galaxies, we highlight the NAT 
galaxy 1330+507 (ID 105). Head-tail and NAT radio galaxies 
are most common in clusters with perturbed X-ray morphologies 
and are probably produced while radio galaxies move through 
cluster atmospheres at high velocities (Burns et al. 1994i. Here 
the NAT galaxy was rejected in the second step of our member 
selection procedure (the galaxy represented by the square with 
the lowest velocity at R~0.7 /i^,,' Mpc in Fig.|5]- top panel) and 
therefore it is likely to be a galaxy "at the border" of the phase- 
space distribution of A1758N. The NAT, projected onto close 



BCGN(SE), points toward A1758N(NW) (see Fig. [B, which 
suggests that it is infalling onto the cluster center. 

In their analysis of galaxies with Spitzer/MlPS 24/im emis- 
sion in the region of A1758, Haines et al. (120091 1 stressed that 
these dusty star-forming galaxies trace the cluster potential and 
in particular the core of A1758N. This supports the hypothesis 
that activity in clusters is triggered by the cluster-cluster merger 
(i.e. is due to the passage of gas-rich galaxies through shocked 
regions of the ICM or affected by time-dependent cluster tidal 
fields, e.g., Roettigeret al. 1996, Bekki 1999; Bekki et al. 20T0^. 
However, they used photometric redshifts to establish the A1758 
membership, while we were able to use spectroscopic redshifts 
to asses their result. We have redshifts for twelve galaxies out 
of infrared luminous galaxies plotted in their Fig. 1 (blue cir- 
cles). Out of these twelve galaxies, five are cluster members and 
they indeed lie in the core of A1758N, between the two BCGs, 
thus supporting the conclusions of Haines et al. (2009; see large 
red diamonds in Fig. [T3T l. Moreover, another four luminous in- 
frared galaxies were rejected from the cluster members only at 
the second step of our procedure, i.e. they lie "at the border" 
of the projected phase-space distribution shown in Fig. |5] (top 
panel). Very interestingly, galaxies with emission lines have a 
similar behavior: all 14 ELGs in our spectroscopic catalog are 
connected to the cluster, eight of them belonging to the cluster 
and six at the border (see Fig.|5]l. In particular, six out of the eight 
members are embedded in the cluster core (see Fig. [13). Thus, 
also the ELGs support the hypothesis that activity in A1758N is 
triggered by the cluster-cluster merger. The effect of the cluster 
merger on the galaxy population was also suggested by Durret 
et al. (201 II). who found an excess of very bright galaxies and 

a bump at Mr' 17 in the observed luminosity function with 

respect to a Schechter luminosity function. 

Here we discuss the relative importance of the NW and 
SE subclusters. From the galaxy clump luminosities, Okabe 
& Umetsu (120081) suggested that A1758N(NW) is the primary 
component and AI758N(SE) is the merging substructure, but 
they measured similar velocity dispersions (crysis ~ 600-700 
km s"') from their weak-lensing analysis. This conclusion also 
agrees with the fact that A1758N(SE) shows an offset with 
the respective X-ray peak and A1758N(NW) does not. Our 
2D-DEDICA analysis shows that A1758N(NW) is richer than 
A1758N(SE), which suggests that A1758N(NW) is the main 
structure. On the other hand, the SE peak shows a comparable 
density (see Table [2]i. This supports the idea that A1758N(SE) 
is the surviving, dense core of a merging structure. The determi- 
nation of the velocity dispersion of subclusters is often a dif- 
ficult task. Figure [T4| - lower panel - shows the cry profiles 
using alternatively BCGN(NW) or BCGN(SE) as centers. As 
for A1758N(NW), the value of its cry is progressively increas- 
ing - which agrees with the small core with peculiar low lo- 
cal cry shown by the DS-test - out to ~ 1000 km s '. As for 
A1758N(SE), the value of its cry is mostly stable at ~ 700- 
800 km s ' before increasing at about half the distance from 
A1758N(NW). Considering the large errors, there is no sig- 
nificant difference between the values of cry of A1758N(NW) 
and A1758N(SE). However, in agreement with the idea that 
A1758N(NW) is the main system, below we consider the nom- 
inal values at the half distance, i.e. cry nw ~ 1000 km s"' and 
o-y,SE ~ 800 km s-i for A1758N(NW) and A1758N(SE), re- 
spectively. 

Making the usual assumptions for the two subclusters (clus- 
ter sphericity, dynamical equilibrium, coincidence of the galaxy 
and mass distributions), we can compute virial global quanti- 
ties. We followed the prescriptions of Girardi & Mezzetti (1200 1 1 
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Fig. 14. Upper and lower panels: integral profiles of mean ve- 
locity and LOS velocity dispersion, respectively. Black and red 
lines consider as centers BCGN(NW) and BCG(SE), respec- 
tively. In both panels the vertical solid lines indicate the distance 
of the BCG of the close subcluster and the dashed lines indi- 
cate the half distance. Thin and thick lines give the results for 
all galaxies and the main system HTll. In the upper panel, the 
two semi-circles give the velocity of the two BCGs. In the lower 
panel, the horizontal line represents the X-ray temperature as in 
Fig.S 



see also Girardi et al. '1998") and computed R^„ - an estimate 
for /?2{)0 - and the mass contained within this radius. In par- 
ticular, we assume for the radius of the quasi-virialized region 
Rvu = 0.17 X o-y/H(z) /i^o Mpc (see Eq. 1 of Girardi & 
Mezz etti |200T] with the scaling of H(z) of Eq. 8 of Carlberg 
et al. 1 1 9971 for /?20())- For the mass we used the virial mass 
corrected for the surface pressure term correction M = My - 
SPT = 3n/2 ■ crlRpy/G - S PT (Eq. 3 of Girardi & Mezzetti 
I2OOII 1. with the size estimate Rpy computed using the full pro- 
cedure based on an assumed typical galaxy distribution and 
SPT = 0.2 X My. In practice, both /?vii and M were computed 
on the basis of the estimated velocity dispersion with the usual 
scaling-laws, where /?vir '^ cry and M(< /?vir) °^ o-y- ^^ com- 
puted Mnw(< Rmr - 2.1 hj^Mpc) - 1.4 xlO'^ /z^Q Mq and 
Mse(< Rvir = 1-7 /z^o'Mpc) = 0.7 xlO'** /i^,} Mq for the NW and 
SE subsystems, respectively. We computed a mass Msys(< R - 
2.1 /i:^qMpc)> 2.1 xlO'^ /z^q Mq for the whole system, depend- 
ing on how far we extrapolated the mass of A1758N(SE) outside 
its Ryir. This value agrees with the value we obtained assuming 
the whole system to be relaxed M(< 7?vii - 2.8 /i^iMpc) = 3 + 1 



Xl0>5 h-^ 



Mp, when rescaled to the smaller R - 2.\ h^ 



Mpc radius. 

As for the cluster structure, finally we discuss our result 
of faint galaxies tracing a peak intermediate between the two 
BCGs, somewhat close to the secondary X-ray peak, while lu- 
minous galaxies trace two subclusters (see Fig. |8). This phe- 
nomenon of luminosity segregation can be connected with the 
cluster internal dynamics. Indeed, very appealingly, galaxies of 



different luminosity could trace the dynamics of cluster merg- 
ers in a different way. A first evidence was provided by Biviano 
et al. ( |I996l l: they found that the two central dominant galaxies 
of the Coma cluster are surrounded by luminous galaxies, ac- 
companied by the two main X-ray peaks, while the distribution 
of faint galaxies tends to form a structure not centered with one 
of the two dominant galaxies, but is instead coincident with a 
secondary peak detected in X-ray data. The observational sce- 
nario of Abell 209 has some peculiar luminosity segregation, 
too (Mercurio et al. 12003) . Therefore, following Biviano et al. 
(1 19961 1, we speculate that the merging is in a post-merging phase, 
where faint galaxies trace the forming structure of the cluster in 
formation, while more luminous galaxies still trace the remnant 
of the core-halo structure of pre-merging subclusters. 



4.2. Merger kynematics and diffuse radio sources 

Although member galaxies can give some pieces of evidence in 
favor of a post-merger phase, the stringent proof is given by the 
X-ray analysis of David & Kempner (2004, Chandra and XMM- 
Newton data) who detected shock-heated gas around the cores 
of the two subclusters. In this study we obtained some impor- 
tant properties of the merger: the LOS relative velocity is $ 300 
km s"' and the merger is a major merger with a mass ratio of ~ 
(2:1). Moreover, in agreement with the results of Pinkney et al. 
(1996), the fact that we clearly detected the subclusters in 2D, 
but not in the ID and 3D analyses, suggests that the plane of the 
cluster merger is mostly perpendicular to the LOS. This would 
also explain the good observability of the X-ray morphological 
features. 

Figure [T] shows the comparison between the positions of the 
density peaks in the galaxy distribution obtained in this work 
and the results from other wavelengths. In agreement with the 
fact that elongated radio halos normally follow the merging di- 
rection (e.g., the "bullet" cluster by Markevitch et al. 2002 , Abell 
520 by Girardi et al. 120081 Abell 754 by Macaiio et al. i20Tll . 
the halo of A1758N is clearly elongated along the NW-SE di- 
rection. As for the two relics, the eastern one (RI in Fig. [T] is 
not perpendicular to the merging axis, as it is typically found 
in double relic clusters (e.g., Abell 3667 Roettiger et al. 1999J 
Kempner & Sarazin 1200 11 1. This supports the idea that A1758N 
is not a head-on merger (see below, but see Bonafede et al. 
I2009l and Boschin et aI. I2010l for alternative explanations of non- 
symmetric radio relics in Abell 2345). More quantitatively, we 
can use our results to examine A1758N with respect to the ob- 
served scaling relation between Pi.4ghz, the halo radio power, 
and the total cluster mass (as determined within ih^Q Mpc, 
Govoni et al. 1200 11 1. Considering the value of f i.4ghz = 9.3 10^"^ 
W Hz ' by Giovannini et al. ( I2009] l and our Msys = 2-3 
xlO'^ /I7Q Mq, A1758N results to be an under-radio-luminous 
cluster (see Fig. 18 of Govoni et al. 2001). However, it should 
be said that A1758N is also very peculiar because of the pres- 
ence of the two close relics. When considering the total diffuse 
emission radio-i-relics Pi.4ghz = 4 10^'* W Hz"' (Giovannini et 
al.|2009j), there is no longer a discrepancy. 

The presence of the diffuse radio emissions suggests that 
A1758N probably is a recent merger because they are expected 
to have a short life time, i.e. on the orderof a few fractions of Gyr 
(e.g., Giovannini & Feretti l2002l Skillman et al. 1201 II I. For radio 
halos, present estimates of the time elapsed after the core-core 
crossing are very short (e.g., 0.1-0.2 Gyr in the "bullet" cluster 
by Markevitch et al. l2U02l 0.2-0.3 in Abell 520 by Girardi et 
al. i2008a . Times are somewhat longer for the radio relics (1 Gyr 
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in Abell 3667 by Roettiger et al. II 9991 1. but they refer to clus- 
ters whose relics are separated by ^ 2 h^^ Mpc, therefore we 
expect that the times for A1758N relics are shorter. An indepen- 
dent evidence of the early age of the cluster merger comes from 
the presence of the bump of highly star-forming galaxies in the 
cluster core, being the duration of the strong starburst induced 
by the cluster merger of < 0.5 Gyr (Bekki 1999). 

On the basis of our results we investigated the relative dy- 
namics of A1758N and A1758S using different analytic ap- 
proaches, that are based on an energy integral formalism in the 
framework of locally flat spacetime and Newtonian gravity (e.g.. 
Beers et al. |1982| l. The values of the relevant observable quanti- 
ties for the two-clumps system are the relative LOS velocity in 
the rest frame, Vi- - 300 km s ' (which is an upper limit); the 
projected linear distance between the two clumps, D ~ 0.8 /i^q 
Mpc; the mass of the system Mjys - 2-3 xlO'^ /z^q Mq. First, 
we considered the Newtonian criterion for gravitational binding 
stated in terms of the observables as V^D < IGM^y^ sin" a cos a, 
where a is the projection angle between the plane of the sky and 
the line connecting the centers of two clumps. The dashed curve 
in Fig. [15] separates the bound and unbound regions according 
to the Newtonian criterion (above and below the curve, respec- 
tively). The system is bound between ~ 3° and 90°; the cor- 
responding probability, computed considering the solid angles 

(i.e., J cos a da), is 95%. We also considered the implemented 

criterion V}D < 2GMsin^ ofy cos a-c, which introduces different 
angles for projection of distance and velocity, not assuming a 
strictly radial motion between the clumps (Hughes et al. |1995l l. 
We obtained a binding probability of 94%. 

In the simple case of a strictly radial motion we can check 
if the merger kinematics agrees with the idea of a very recent 
merger because the simple linear-orbit two-body model connects 
the time elapsed from the the core-core passage t and the pro- 
jection angle a. Figure [15] shows the bimodal-model solutions 
as a function of a. The thick/black lines show the case with 
Vrf.LOS = 300 km s^'and the time from the core-core passage 
t - 0.2 Gyr. There is a bound outgoing solution (BO) with a ~ 
15-20° that leads to a deprojected relative velocity of Vrf ~ 
1150-900 km s"'. The effect of decreasing the relative veloc- 
ity is that of decreasing the projection angle, e.g. Vrf.LOS - 100 
km s^' leads to a ~ 5-10° and thus to Vii> 1 150-600 km s^'. 
The effect of increasing the age of the merger is that of increas- 
ing the projection angle and decreasing the relative velocity, 
e.g. t - 0.5 Gyr leads to a ~ 55-60° and thus to Vd ~ 350 
km s"' (here there are also two symmetric bound incoming so- 
lutions). Therefore only a very recent merger agrees with the 
evidence that the merger plane is mostly perpendicular to the 
LOS. Moreover, we have to consider the limits on the relative 
velocity indirectly recovered from X-ray data. In order to pro- 
duce the shock waves, which e.g. originate the relics, we need a 
Mach number M > I and David & Kempner ( 20041) computed 
A1$ 1.15 from X-ray data. The Mach number is defined to be 
M - Vs/cs, where Vj is the velocity of the shock and c^ in the 
sound speed in the pre-shock gas (see e.g., Sarazin 120021 for a 
review). Assuming for the sound speed in the pre-shock gas the 
value obtained from the X-ray temperature of David & Kempner 
(6-9 keV), i.e. v^ = 1000-1200 km s"', and v^ ~ Vrf, we expect 
Vrf to be in the range of 1000-1300 km s ', which agrees with 
models with a very recent merger. 

The obvious limit of the model described above is that the 
merger is treated as a head-on merger. Although general argu- 
ments suggest small impact parameters for cluster mergers ( ~0.2 
hj^ Mpc, see Sarazin |2002| |. it seems that this is not the case for 
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Fig. 15. System mass vs. projection angle for bound and un- 
bound solutions (thick solid and thick dashed curves, respec- 
tively) of the two-body model apphed to the A1758N(NW) 
A1758N(SE) subclusters. Thick/black Unes refer to t = 0.2 Gyr 
and Vrf.LOS - 300 km s ' . Labels Bla and Bib indicate the bound 
and incoming, i.e., collapsing solutions (solid curve). Labels BO 
and UO indicate the bound outgoing, i.e., expanding solutions 
and unbound outgoing solutions (solid curve going on in the dot- 
ted curve, respectively). The horizontal lines give the range of 
observational values of the mass system. The dashed curve sep- 
arates bound and unbound regions according to the Newtonian 
criterion (above and below the thin dashed curve, respectively). 
The red/faint lines refer to f = 0.5 Gyr and Vrf.LOS = 300 km s"' . 
The green/faint lines refer to f = 0.2 Gyr and Vrf.LOS = 100 
km s"'. 



A1758N. From the morphological details of the X-ray surface 
brightness image, David & Kempner ( I2004I I deduced that the 
NW subcluster is moving north and the SE subcluster is moving 
southeast and that A1758N is in the late stage of a large im- 
pact parameter merger: their Fig. 9 shows the proposed orbits 
in the merging plane. However, these orbits are not supported 
by the two elongated regions of high metallicity likely left by 
ram-pressure stripping during the merger (Sect. 4 and Fig. 1 3 of 
Durret et al. '201D: the NW feature is elongated along NW-SE 
as in the case of a head-on merger, while the SW feature has a 
clear tail perpendicular to the NW-SE direction as in the case 
of an off-axis merger. Thus, to date, we are far from knowing 
the details of the possible orbits and from investigating a more 
complex kinematics. 

Our general conclusions about A1758 agree with other typ- 
ical clusters that show extended diffuse radio emissions, i.e. we 
are dealing with a massive cluster with a major ongoing merger 
(like, e.g., Abell 2744; see Boschin et al. 120061 ) and, indeed, 
less massive merging systems do not seem to host extended ra- 
dio emissions (e.g. Abell 2146 by Russell et al. 120111 see also 
Rodriguez-Gonzalvez et al. 1201 II ). 
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